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In the spectator model of heavy hadron decay, the Ufe- 
times of aU 5-hadrons are equal, independent of the fla- 
vor of the lighter quarks bound to the b quark. Using 
the heavy-quark expansion [H 01 in the calculation of 
the width, spectator quark interactions enter in higher 
order {AQco/'mb)'^ terms where mf, is the mass of the 
b quark and Aqcd is the energy scale of the QCD in- 
teractions within the hadron. This leads to the lifetime 
hierarchy t{B'^) = t{B^) < t{B^). Theoretical results 
predict t{B+)/t(B°) = 1.06 ± 0.02 and t{B°)/t{B°) = 
1.00 ± 0.01 0, Si- The world averages for the corre- 
sponding experimental numbers are 1.071 ± 0.009 and 
0.965±0.017, respectively Q. The precision of our knowl- 
edge of the i3° lifetime is much less than for the B^ and 
B~^ lifetimes, and therefore, a more precise measurement 
would be useful, both in general and for comparison with 
theoretical calculations. Such a measurement is espe- 
cially warranted since the agreement on the lifetime ratio 
between theory and experiment is only fair. 

In this Letter, we present a measurement of the life- 
time in flavor-specific decay modes. The data come from 
pp collisions at ^/s = 1.96 TeV at the Fermilab Tevatron. 
This analysis is based on an integrated luminosity of 1.3 
fb~^ collected by the CDF II detector between Febru- 
ary 2002 and November 2006. This sample yields more 
than 1100 fully-reconstructed i?j — D~n^ candidates 
with Dj — >■ (j)TT~ and (p — )■ K^K~ after online and of- 
fline selection @. In addition, the sample reconstructed 
as — )■ DjTr~^ includes partially-reconstructed (PR) -B^ 
candidates that are used in this lifetime measurement and 
more than double the number of 5° candidates available 
for analysis. One such PR decay is -B° — >■ D~ p'^ with 
p"*" — >■ tt+tt" where the tt^ is not reconstructed. The in- 
clusion of PR decays introduces an uncertainty in the 
momentum measurement of a given candidate. However, 
a correction to the proper decay time has been estimated, 
and the total uncertainty on the lifetime measurement is 
improved by the use of the PR final states. 

The CDF II detector is described in detail in Ref. 0- 
The detector elements relevant for this analysis are the 
silicon vertex detectors and the central drift cham- 

ber (COT) [llj. The silicon detectors consist of 7 or 8 
layers of microstrip silicon sensors covering the pseudo- 
rapidity range \r]\ < 2.0. The COT is an open cell 
drift chamber covering {rjl < 1.0. Both the COT and 
silicon vertex detectors are immersed in a uniform 1.4 
T axial magnetic field with the field axis parallel to the 
proton beam. 

A data sample enriched in hadronic B decays is se- 
lected with a three-level trigger system that searches for 
tracks displaced from the primary vertex. At level 1, pat- 
terns of hits in the COT are identified as tracks by the 
extremely fast tracker (XFT) [13]. At level 2, the sili- 
con vertex trigger [l^l associates a set of silicon hits with 
the XFT tracks and improves track measurement preci- 
sion. The trigger requires each event to contain a pair of 
charged particle tracks, each having transverse momen- 
tum pt ^ 2 GeV/c and transverse impact parameter do 



in the range do € [120 fim, 1 mm], where do is defined 
as the distance of closest approach between the particle 
trajectory and the beamline, measured in the transverse 
plane. The opening angle between the tracks' trajecto- 
ries (A0 in the plane transverse to the beam) must be be- 
tween 2° and 90°, and their intersection must be at least 
200 /im from the interaction point, as measured in the 
plane transverse to the beam direction. At level 3 track 
reconstruction is performed entirely in software, with the 
full precision of the tracking system available, and the 
level 1 and 2 requirements are confirmed. These trigger 
requirements preferentially select events containing long- 
lived particles and sculpt the proper time distribution of 
the particles that are accepted for analysis. CDF also em- 
ploys two more restrictive triggers that require the tracks 
in the trigger pair to have opposite charges, individual 
Pt > 2(2.5) GeV/c, and the scalar sum pr > 5.5(6.5) 
GeV/c. 

We reconstruct _B" — D'^tt^ candidates (where B^ 
and 1^7 imply i?° candidates and candidates) 
by first identifying D~ — >■ (j){K^ K^)tt^ from tracks 
with Pt > 350 MeV/c using the invariant mass re- 
quirements \m{K-K+) - 1020. 5| < 7.5 MeV/c^ and 
\m{K-K+T:~) - 1968.3| < 20 MeV/c^. The Dj daugh- 
ter tracks must satisfy a three-dimensional vertex fit. We 
then combine each Dj with a positively charged track 
with Pt > 1.0 GeV/c to form a iJ^ — )■ D~tt^ candi- 
date and require the pair to satisfy an additional three- 
dimensional vertex fit. We do not constrain the mass of 
the (j) or Dj in this fit. The decay length of the is 
measured with respect to the event's primary vertex and 
must satisfy requirements on the following quantities: the 
decay length of the 5° projected along the transverse 
momentum, Lxy{B^) > 450 /im, and its significance, 
Lxy{Bs)/o'L^y{Bg) > 5; the transverse distance between 
the B^ and D~ decay points is greater than 0; the trans- 
verse impact parameter of the 5°, |do(-B°)| < 60 /im; 
and the significance of the longitudinal impact param- 
eter, \zoiB°)/a^g{Bf)\ < 3. Both fits for the 5° and 
D~ vertices must have reasonable goodness-of-fit values 
when considering only the track parameters measured in 
the transverse plane. 

To further separate B^ mesons from backgrounds with 
similar topologies, we require the transverse momen- 
tum of the Pt{Bs) > 5.5 GeV/c, and the an- 
gular separation between the Dj and the tt from the 
fiO, AR{Dj,ttb) = a/(A77)2 + (A0)2 < 1.5. We re- 
quire the isolation of the 5° to be greater than 0.5, de- 
fined as pt{B1) divided by the scalar sum of the trans- 
verse momenta of all the tracks in a cone of Ai? < 1 
around the B^. We tighten the requirement on the 
mass of the Dj (|m- 1968.3| < 12 MeV/c^) and veto 
Dj candidates consistent with D*^ — > D^tt^ where the 
decays to K^n^ii^ . The D* veto is accomplished 
by taking the D~ daughter tracks {K~ K'^'k~), assign- 
ing the pion mass to the negative kaon, and requiring 
Am = m{K+'!T~TT') - ra{K+TT-) > 180 MeV/c^. We 
also require that the decay contain two reconstructed 
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tracks satisfying the level 2 trigger requirements. 

The simulated data samples used in this analysis con- 
sist of single 6-hadrons generated by BGENERATOR [H, 
[l6j with pt spectra consistent with NLO QCD and de- 
cayed with EVTGEN jTTi]. Full detector and trigger sim- 
ulations are performed. The simulated B candidates are 
reconstructed with the same procedure and the same se- 
lection as the data candidates. We reweight the simu- 
lated sample to match the data distributions for pt{B) 
and trigger scenario mixture. 

The lifetime of the meson is determined from two 
sequential fits. The first is a fit to the invariant mass 
distribution of candidates reconstructed as Djn'^ used 
to determine the fractions of the total number of events 
found in the various decay modes. These fractions are 
then used as fixed inputs in the second fit to the proper 
decay time distributions of the B^. The uncertainties 
on the fractions returned by the mass fit are treated as 
sources of systematic uncertainty. 

The mass fit is an unbinned maximum likelihood fit 
to the invariant mass of the candidate reconstructed as 
Djn+ with m^'j^'' € [4.85,6.45] GeV/c^ The mass fit 
components can be characterized as coming from one of 
three possible sources: single &-hadrons, real-I?J-f track 
background, and fake-Dj-|-track background. The mass 
probability distribution functions (PDFs) for single b- 
hadrons were obtained from simulation, with an addi- 
tional small shift and resolution smearing to bring the 
simulated -B^-peak central value and width into agree- 
ment with data. The single-^ modes were separated in 
the fit as follows: ^ D-TT+{n-f), B2 D^K^, 



B" 



D: 



D7*n+. B° -> Ds ^*^X, B°/B- 



D-X, B"^ ^ D 



TT' +Ds^*'^K+, and A,, ^ A^X. 
ratio was constrained to the results 



The DfK^/D- 
of Ref. [18]. 

Real-Dj-l-track backgrounds consist of a real I?^, pro- 
duced promptly or from a &-hadron decay, plus an ad- 
ditional track produced in the event. The mass PDF 
for these events is obtained from an auxiliary fit to the 
wrong-sign sample, which consists of data events recon- 
structed as D~'K~ and sideband subtracted in the Z?^ 
mass. The mass PDF for the fake-D^" -|-track background 
is obtained from an auxiliary fit to the Dj sidebands. 
The results of the mass fit are shown in Fig. [1] with 
various modes combined for plotting only. The real- 
Dj-l-track background and fake-ZJj-l-track background 
components are drawn together as the "combinatorial 
background" . 

For the lifetime fit, the variable of interest is the proper 
decay time, defined as ct = {L^y{B^^)-m)^'')/pT{B^,). The 
reconstructed mass vnT^'^ is used instead of the world av- 
erage B^ mass. A salient feature of this analysis is the 
treatment of partially-reconstructed mesons as sig- 
nal events that contribute to the lifetime measurement. 
Since in the partially-reconstructed cases Lxy{B^), m^g'^, 
and pxiB'^) are extracted from candidates that are miss- 
ing particles after reconstruction or have the wrong mass 
assignment for a daughter particle, a multiplicative cor- 
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1: Mass distribution for candidates reconstructed as 
Dj-K^ with fit projections overlaid. 



rection factor K to the decay time is needed. K is defined 



K 



[pt(S°) 



J 



i^'^'^-cos^ypRj 



where 9pn. is the angle in the x—y plane between the true 
momentum of the B^ and the momentum of the partially- 
reconstructed Because the ratio nf^'^/pT{B^) is nu- 
merically very close to the ratio mg'"'^/p^"^(B°), this 
choice of ct definition forces the K factor distributions 
to be centered near K = 1 with widths of a few percent. 
The K factor distributions are determined with simula- 
tion. 

The lifetime of the B^ meson is determined from an un- 
binned likelihood fit to the B^^ candidates with invariant 
masses in the range [5.00, 5.45] GeV/c^. There are three 
main types of lifetime fit components that will be de- 
scribed in the following paragraphs: fully-reconstructed 
Sg, partially-reconstructed 5°, and backgrounds. The 
treatment of each component depends on its decay struc- 
ture and whether it can provide information about the 
S° lifetime. 

Fully-reconstructed (FR) modes where all of the B^ 
daughter particles are included with the correct mass as- 
signment in the construction of the B^ candidate are the 
first type of lifetime fit component. The only FR mode 
in this analysis is the DJtt^. The core functional form 
of the FR PDF is an exponential with decay constant 
ct{B1) convoluted with a Gaussian resolution function 
with width a: 



1 



-ct' 



1 



2'Ka 



eS{ct) (1) 



A multiplicative "efficiency curve" accounts for the trig- 
ger and analysis selection criteria: 
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eff(rf) =^iV^•(c^-/?^)' 



Fig. [2l The results of the fits performed separately in the 
FR mass region (1.456 ± 0.067 ps) and PR mass region 
(1.544 ± 0.051 ps) agree with each other at a level of 1.0 



The shape parameters (cr, /3i, Ni and r^) of the PDF 
are determined in a fit to a simulated i?^ sample where 
the lifetime used for generation is known. All the pa- 
rameters for the PDF are then fixed and only T(i?°) is 
varied in the final fit to the data. As we depend on the 
simulation of the displaced-track trigger, we use a data 
sample of J/?/' — > M^A*^ decays collected with a di-muon 
trigger to assess the accuracy of this assumption and as- 
sign a "trigger simulation" systematic uncertainty based 
on these studies. The partially-reconstructed, PHOTOS- 
modeled D~Tr{n'-f) decays [l^ are combined with the FR 
DjTT, as the momentum carried by the photon is small 
and their lifetime distribution is extremely close to the 
FR one. This simplification is considered as a possible 
source of systematic uncertainty. 

Partially-reconstructed modes either neglect J5° 
daughter particles in the construction of the B'^ candi- 
date or assign them an incorrect mass. B'^ — > D^K^, 



D: 



d: 



and other decay modes partially recon- 



structed under the Djn'^ hypothesis can also contribute 
to the i?5 lifetime measurement. The PR PDF is similar 
to the FR PDF of Eq. ([T]) with an additional convolution 
with the K factor distribution for each mode. There are 
separate efficiency curve parameters for each mode, again 
determined from fits to simulated events. 

The backgrounds in the lifetime fit can either come 
from decays of 6-hadrons other than the 5° {e.g., 
B°/B- D-X, 5° ^ DjX, and A,, -> AcX), or 
they can come from real-_D^+track and fake-Z?^ +track 
combinations. The PDFs for the former modes are de- 
rived from fits to simulated i?*^, B~ and A;, samples, 
while the ones for the latter combinatorial backgrounds 
come from the two proxies available: the -B^ upper side- 
band taken from the m^^'^ interval [5.7, 6.4] GeV/c^ and 



sidebands. The sidebands are taken from 
interval [5.0, 5.45] GeV/c^ and the to^'^ in- 



the D; 
the m\ 

terval [1.924,1.939]U[1-999,2.014] GeV/c^. Both prox- 
ies contain a mixture of fake Dj-l-track events and real 
Dj-|-track events, where a real Dj can be poorly recon- 
structed. All the background shape parameters are fixed 
in the final lifetime fit. 

The analysis procedure was tested extensively on three 
control samples: B'^ — > D^tt^ with — > K^tt^tt^, 
B" -> D*-7r+ with £)*-_-> D"tt- and D° K+n-, 
and 5+ — £)°7r+ with D° — > A'+tt^ before perform- 
ing the Bg fits. Furthermore, the lifetime fit of the B^ 
was performed using a blind approach, i.e. by determin- 
ing the statistical and systematic uncertainties without 
knowledge of the fit result itself. Good agreement with 
the world average values [Bl of the B^ and lifetimes 
was found. 

The lifetime of r(B^) = 1.518 ± 0.041 (stat.) ps is 
obtained from the full fit. The fit results are plotted in 
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FIG. 2: Distribution of ct for candidates reconstructed as 
Bg — >■ DJ {(fnT^)Tv'^ with fit projection overlaid. 



TABLE I: Summary of sources of systematic uncertainty for 
the Bg — > DJ {4>-K^)X lifetime fit. The total uncertainty is 
calculated assuming the individual contributions are uncorre- 
lated. 



Description 



Value (ps) 



Background modeling and fractions 0.019 

Fixed single-6 background ct 0.003 

Reweighting for pT and trigger 0.012 
Lifetime contribution of Dtt radiative tail 0.002 

Efficiency curve parameterization 0.002 

Trigger simulation 0.014 

Impact parameter correlation 0.003 

Detector alignment 0.003 



Total systematic uncertainty 



0.027 



We use a Monte Carlo technique to assess the system- 
atic uncertainties. For each source of systematic uncer- 
tainty, we generate 1000 simulated experiments with the 
number of events in each experiment Poisson-distributed 
around the number of events in data. The simulated ex- 
periments are generated with a non-standard lifetime fit 
configuration (where the PDFs or numbers of events in 
the various modes are modified to account for the system- 
atic effect) and fit with the default configuration. The 
mean biases returned from the fits to the simulated ex- 
periments (Tret — Tgcn) are used to set the size of the 
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systematic uncertainties. We consider several sources of 
systematic uncertainty: combinatorial background frac- 
tion, modeling of backgrounds from single 6-liadron de- 
cays, effect of reweighting the full simulations to match 
the data, modeling of the trigger bias as a function of 
ct, offline-online impact parameter correlation, accuracy 
of the trigger simulation in Monte Carlo, and detector 
alignment. Table |T] contains the final list of systematic 
uncertainties for this measurement. The largest contri- 
bution comes from the uncertainty on the total amount of 
combinatorial background and the amount of promptly 
produced real-Z?^ background. 

The displaced-track trigger, in addition to modifying 
the accepted decay length distribution from a simple ex- 
ponential to the form in Eq. ([1]), alters the expected 
mixture of mass eigenstates i?^ ^ and B'^ ^ in the fiavor- 
specific i?" — > D^TT^ decay by preferentially selecting the 
longer lived i?" . The size of the imbalance can be cal- 
culated using the parameters of the efficiency curve and 
the world average of l/F = 1.47 ps Our result can be 
corrected back to a flavor-specific lifetime measurement 
with (5t(B°) = -0.11 {AT/rf ps. Given the world av- 
erage Ar/r = 0.092j:°;°^^ [|, the correction would be 
smaller than our statistical and systematic uncertainties. 
Therefore, we do not correct the central value or assess 
an additional systematic uncertainty. 



In summary we have measured the B'^ lifetime using 
both fully reconstructed i?° — )• -D7(07r~)7r+ and partially 
reconstructed B'^ _D7(07r~)X decay modes, in a sam- 
ple with 1.3 fb^^ of integrated luminosity. We measure 
t(BO) = 1.518 ± 0.041 (stat.) ± 0.027 (syst.Vps, which is 
consistent with theoretical expectations 0, IJ- 
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